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Abstract 

Given a joint probability density function of N real random variables, obtained 

from the eigenvector-eigenvalue decomposition of N x N random matrices, one constructs a 
random variable, the linear statistics, defined by the sum of smooth functions evaluated at the 
eigenvalues or singular values of the random matrix, namely, 

For the jpdfs obtained from the Gaussian and Laguerre ensembles, we compute, in this 
paper the moment generating function E^(exp(—A Yj where denotes expectation 

value over the Orthogonal (/3 = 1) and Symplectic (/3 = 4) ensembles, in the form one plus a 
Schwartz function, none vanishing over M for the Gaussian ensembles and M'*' for the Laguerre 
ensembles. 

These are ultimately expressed in the form of the determinants of identity plus a scalar 
operator, from which we obtained the large N asymptotic of the linear statistics from suitably 
scaled F{-). 


1 Introduction 

The well-known joint probability density function for the eigenvalues of N x N Hermitian 

matrices from an orthogonal ensemble (d = 1), unitary ensemble {(3 = 2) or symplectic ensemble 


1 


(/9 = 4) is given by [H] 


N 


N 


P^^\xi,X2,... ,XN)YldXj = JJ 


Xj — Xkf w{xj)dxj, 
j=l l<j<k<N j=l 


where w(x) is a weight function or a probability density supported on [a,b]. In this paper, we 
require that the moments of w, namely, 


Sj := / x^ w{x)dx, 


to exist for j = 0,1,2,.... Here the normalization constant reads, 

4/3) 


r<\p) _ 


/[a,b]iv ni<j<fc<V \^j - ^kf nf=l w{Xj)dXj 
Many years ago, Selberg |T2], obtained closed form expression for where 
w{x) = x“(l — x)^, a > —1, b > —1, X G [0,1], the Jacobi weight. The constant with the 
Gaussian weight w{x) = , x G M and Laguerre weight, w{x) = x" e”"^, a > — 1, x G can 

be found in HI]. 

If we take w{x) = , x G M, these are known as the Gaussian orthogonal ensemble (GOE), 

Gaussian unitary ensemble (GUE) and Gaussian symplectic ensemble (GSE). If 
w{x) = x“e“^, a > —1, X G M+, then we have analogously the LOE, LUE and LSE. 

The moment generating function, 

g(J)(/) - Efl (e-' ■= nia<t<A' - “^‘1' nf.i toixj) e-" 

is given by 


ni<j<fc<V l^i ^k\ Y[j=l''^i^j)dXj 


N 


G'!v^(/) = / n \xj - Xkff\w{xj)[l +f{xj)]dxj, (1.1) 

i<,<fc<jv i=l 

where /3 = 1,2,4. Here ;= 1 + /(x), and we do not indicate that /(■) also depends on A. 

Throughout this paper, we assume /(x) lies in the Schwartz space, and 1 + /(x) ^ 0 over [a, b]. 

The simplest, very well-studied, unitary case corresponds to j3 = 2, see [D HI El [31 El E], and 
the references therein. 

We state here, for later development, facts on orthogonal polynomials. Let 
(Pj{x) := Pj{x)^yw{x), j = 0,1,2,..., where Pj{x) are the orthonormal polynomials of degree j 
with respect to the weight w{x) supported on [a,b], 


Pj{x)Pk{x)w{x)dx = Sj^k, j,k = 0,1,2,.... 


2 





It is a well-known fact that, 


G®(/) = det (/ + A'®/) . 

where f is an integral operator with kernel K^^\x, y)f{y)- Here K^^\x, y) := Yl!j=o 
can be evaluated via the Christoffel-Darboux formula. 

With Gaussian background, where w{x) = e“^^, x G M, we have 

, c, =7r32VrO' + l), J = 0,l,2,..., 


PM) ■■ = 


where Hj{x) are the Hermite polynomials of degree j. 

With Laguerre background, where w{x) = x°‘e~^, a > — 1, a; G M’*', we have 


- (a) 


R{x) := 


M (a) _ /r(j -|- Q -|- 1) . _ 




= 


r(j +1) 


, J = 0,l,2,..., 


where L^^\x) are the Laguerre polynomials of degree j. Properties of the Hermite and Laguerre 
polynomials can be found in [T3] . 

We shall mainly deal with the (3 = 4 and (3 = 1 cases, which are more complicated than (3 = 2 
situation, especially so for /3 = 1. 

We always begin with the general case, and then apply the results obtained to the two classical 
ensembles, i.e., the Gaussian ensembles and Laguerre ensembles. 

Furthermore, in (3 = 1 situation, for convenience, N is taken to be even, and it is expeditious 
to make use of the square root of the Gaussian weight, x G M and the square root of the 

Laguerre weight, x°‘Pe~^P^ a > —2, x G M"*” in later discussion. 

In this paper, we shall be concerned with the large N behavior of for the Gaussian 

ensembles and Laguerre ensembles. For the Gaussian ensembles, we replace f{x) by / {^2Nx^ , 
in the orthogonal case, and / , in the symplectic case. For the Laguerre ensembles, we 

replace f{x) by / ^\/4A^xj, in the orthogonal case and / ^\/8A^xj, in the symplectic case. 

We will ultimately give the mean and variance of the linear statistics, as N ^ oo, together with 
leading correction terms. 

For comparison purposes, we write down results in the GUE, where, w{x) = , x G M. 

Denote by and the mean and variance of the linear statistics P ^\/2iVxjj. It 

is shown in [3], that, 

1 

F{x)dx, N oo, (1.2) 


[GUE) 

EN 


^ 1 
TT 
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V 


(GUE) 

N 


1 

-)■ - 
TT 


^oo /*oo 


F‘^{x)dx 


sin(a; — y) 


n 2 


F{x)F(y)dxdy, N ^ oo. (1.3) 


—ooJ—oc _ y) 

For the LUE, where w{x) = x‘^e~^, a > — 1, x G M+. Denote by and the 

and variance of the linear statistics ^ (-^diVxj). It is shown in[2], that, 


mean 


poo 

^{LUE,a) ^ / B‘^'^\x,x)F{x)dx, iV —t oo, 

Jo 


■,XLUE,a) 

where 




B^^\x, x)F‘^{x)dx 


(x, y)B^^^ (y, x)F{x)F{y)dxdy, 


B^^\x,y) : = 
B^°‘\x,x) : 


JMyJai.y) - Jc{,y)xJ'^{x) 
- 2 - 2 -^ 

Jl{x) - Ja_i(x)J„+i(x) 

= -::- X. 


(1.4) 
t oo, (1.5) 


( 1 . 6 ) 

(1.7) 


Here Ja(-) is the Bessel fnnction of order a. 

We want to point out, the motivation of this paper comes from [B], which provided results both 
for symplectic ensembles and orthogonal ensembles, and specialize to the Gaussian case, i.e., GSE 
and GOE. The [6] dealt with the situation where /(.) is the characteristic function of an interval 
(or the union of disjoint intervals), and focus on the distribution of the mth largest eigenvalue in 
the GSE and GOE, while we are interested for ’’smooth” /, and we also consider the Laguerre case, 
i.e., LSE and LOE. 

This paper is organized as follows. In Section 2, we recall a number of theorems, the operators 
D and £, and end with two Lemmas relevant for later development. Section 3 begins with a general 
discussion of the Symplectic ensembles in a general setting, followed by detailed discussions on the 
GSE and LSE cases and ends with the computation of the mean and variance of linear statistics for 
large N. Section 4 repeats the development in Section 3 but for the Orthogonal ensembles, which 
is harder. We conclude in Section 5. 


2 Preliminaries 

For orientation purposes, we introduce here a number of results, which will be used throughout this 
paper. 
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Theorem 2.1. The Stirling’s formula fT^ 


T{n) = y/^n’^ 2 e ” [l + O (n , n —)■ cxd. 


( 2 . 1 ) 


Lemma 2.2. 


2f. (-/V, 1 ; 2) = + o (Af-‘). Af ^ 

Proof. From the integral representation of the hypergeometric fnnction, 

2^1(0,/?;7; z) = ^ Re7 >Re/ 5>0, 

we obtain, 


2-^1 




\/r^ 


Let 

X = \/l — t, 


( 2 . 2 ) 


then 


2R1 




l)'^(ix 



l)^fix + 



l)^dx 


- 1 ) 


N 


.LI 

2 


(1 - 2x^)^dx + 



l)^(ix. 


(2.3) 


Consider the hrst integral in fl2.3p . Let 


\/2 

X = cos 9, 9 G 


TT' 


2J ’ 


then 


• LI 


(1 — 2x^)^ dx = 




(sin 9) 


\/2 [2'^r(iV + 1)] 


2 


2 r(2Ar + 2) 


where nse has been made of the Stirling’s formnla (12.ip . in the last eqnality. 
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Now consider the second integral in fl2.3p . 


[ {2x^ — 1)^dx < [ {2x — l)^dx 

’ V2 I v/2 

2 ^2 


Hence 




l-(\/2-1)^+1 

2{N+1) ■ 


Lemma 2.3. If a > 0, then 

2n / 




1)’ 


2n + a 


UZJ» + 20 \2n-m + ij {2n + 1)! 
where (1) := -!Nz±l)- 


, n = 0,l,2, 


r(k+i)r{j-k+i) ■ 


Proof. Firstly, 


2n 

E 

m=0 


■1)' 


m 


!n£m(« + 20 \2n-m + l 


2n + a 


1 r(2?7, + a + l)r(l) 2 Ei( - 2?7, - 1, a; 2) 


(2n + 1)! 22"+ir(2n + 2)r(2n + f + l)r(a) 


From fl2.2p . we see that, 


..n ( - 2n - 1, a; 2 ) = - () t-i(l - 


Let 


then 


t = cos^ ( ^ )) ^ ^ [0, tt]. 




OL—1 


(sin0)“-^(cos0)2"+M0 


□ 


(2.4) 


We carry ont mathematical indnction in n to prove that for a > 0, fj(sind)'^ ^(cos6*)^"'+^(i0 = 0. 
Setting if n = 0, in the integral above, we have, 

f (sin0)“-^cos6'd0 


sin6*)“ ^dsin6' 


(sin0)^ 


a 


= 0 . 
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Suppose /J^(sin6*)" ^{cos d9 = 0, for a > 0, then 

rn 

/ (sin0)"-i(cos6»)2^+3(i0 

Jo 

= [ {sme)^-\cosef^^^cos^9de 

Jo 

= / (sin 9)°-^ (cos (1 - sin^ 9) d9 


{sin 9)^-\cos 9Y^^^d9 - I {sin 9)^+\cos d9 


0+1 / 


\2fc+l 


= 0 . 


Hence, 


It follows that, for a > 0, 


(sin0)“-^(cos0)""+M0 = 0, a>0. 


2 Fi( -2n- l,|;a;2) = 0, 


and from fl2.4p . the lemma follows. 


□ 


Theorem 2.4. 

TT {xj - Xk)'^ = det (xl,jxl~^)j^o,..., 2 N-i , 

l<j <k<N 

where the determinant on the right is a 2N x 2N determinant with alternating columns 
{x\}, {jx\~^}, {x^}, {jx^-^},... If. 

The next theorem, due to de Bruijn [5], is instrumental for the hnite N computations in Sections 
3 and 4. 


Theorem 2.5. For any integrable functions Pj{x) and qj{x), j = 1,2, we have 


[a,b]^ 


det {pj{xk),qj{xk))j=h..., 2 N dxi--- dx^ = {N\Y det / {pj{x)qk{x) - Pk{x)qj{x)) dx 


2N 


j,k=l 


det {Pj{xk))fk=i dxi--- dxN ] = det ( 

I / V 


b rb 



N 


sgn{y - x)pj{x)pk{y)dxdy j 

/ j,k=l 


' a<xi<-"<XN'^b 

where the determinant on the left of the first equality is a 2N x 2N determinant with alternating 
columns {pj{xi)}, {^^(xi)}, {pj{x 2 )}, {(lj{x 2 )}, ■ ■ ■■ In addition, N is even in the second equality. 
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Theorem 2.6. If A, B are Hilbert-Schmidt operators on a Hilbert space H, then & 


det{I + AB) = det{I +BA). 

Following [HI El E], we introduce here operators e and D which will be crucial for later 
development. Let e be the integral operator with kernel 

:= ^sgn(a:-?/), 

then for any integrable function g dehned on [a, b], 

= j e{x,t)g{t)dt=^(^j g{t)dt- j g{t)dt^ , xe[a,b]. 

It is clear that s{y,x) = —e{x,y), i.e., £* = —e, where * denotes transpose. 

Let D be the operator that acts by differentiation, thus for any differentiable function g dehned 
on [a, h], 

Dg{x) = = g'{x). 

With further conditions on g{x), we prove an easy lemma on the commutator [D, e\. 

Lemma 2.7. For any function g G C^[a,b] and g{a) = g{b) = 0, Deg{x) = eDg{x) = g{x), i.e., 
De = eD = I. 

Proof. For any function g G C^[a, b] and g{a) = g{b) = 0, we have 

{De)g{x) = -^ f e{x,t)g{t)dt 


2 dx 


g{t)dt — / g{t)dt 


= + \ 9 {.x) = g{x), 


and 


{eD)g{x) = / e{x,t)g'{t)dt 

J a 

1 r 1 

- 2 / 


g'{t)dt 


= \[ 9 ix) - 9 ( 0 -)] - ^[gib) - g{x)] 
= 9{x) - ^[9{a) + 9{b)] 

= 9{x). 





The proof is complete. □ 

Denote by u®v the integral operator with kernel {u®v){x, y) := u{x)v{y). We have the following 
lemma. 

Lemma 2.8. If A is an integral operator with kernel A{x,y), then 

A{u 0 v) = (Au) 0 V, (u0 v)A = u0 (A^v). 


3 The symplectic ensembles 


3.1 General case 


Taking /3 = 4 in fll.ip gives the symplectic ensembles and generating function becomes, 


Gf(f) = C 


(4) 

N 


f 4 ^ 

/ TT {Xj - Xk)^T[w{xj)[l + f{xj)]dxj, 

l<j<k<N j=l 


where 


_ 

— 


Xk) Ilj=i'w{xj)dxj 

is a constant depending on N. 

We follow closely the computations of Dieng [6], and Tracy and Widom [15 
From Theorem 12.41 and Theorem 12.51 and some linear algebra, we get 


G'!v^(/) = det ( / [nj{x)7r'f.{x) - 7r'-{x)7rk{x)] + f{x))dx 


44 ) 


27V-1 


j,k=0 


where is a dependent constant, and TTj{x) is any polynomial of degree j. Let 

ijj{x) = 7lj{x)^/w{x), 

and following PE], we see that, 

GS^^(/)]' = det 

where are matrices given by 


= ( / {'tpj{x)'ijj',^{x) - ^jJj{x)'^jJk{x)) dx 


2N-1 


j,k=0 


(3,1) 


9 










With the notation, 




{ijj{x)^k{x) - ^jJJ{x)^jJk{x)) f{x)dx 


we obtain, hnally. 




W(/)]" = det + 


where is the integral operator 


^(4) _ 


Y.fk=l d'jk'ipj - E?J=o ® 

®i^k - E?S ^jki^j ® fk 


/y(4>4) /yb.2) 

-''-AT,4 -'V^4 

/y(2.1) 7^(2,2) 

AT,4 -'AJY4 


In the next theorem, we obtain further relations on = 1,2, which ultimately expresses 

in terms of . 


Theorem 3.1. 


/y(2.1) _ /y(2:2) /yb,l) _ j^t^{2,2) /y(4,2) _ j^-^(2,2) ^(2,2) p. _ ^(2,2) 

A^,4 A^,4 ' 


Proof. First of all. 


j^na ^ 




j,k=0 

2N-1 

Y d'jkf’j ® i^k 


~ -''■AT,4 • 


Secondly, for any integrable function g{x) dehned on [a,b], we have 


DK^NAg{x) = ^ K^NA{x,y)g{y)dy 


-k^na i.^^y)9{.y)dy 


= / d^NA{^^y)9{y)dy 

J a 

= Kna9{x), 
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i.e., 


It follows that 


Similarly, 


~ -''■AT,4 • 




DK^NAg{x) = ^ I K^Nf{^^y)9{y)dy= I -^K^Ni\^^y)9{y)dy 


d 


Kna {^^y)9{y)dy 


= d<^NA9{x), 


I.e., 


Finally, we find 


^(1-2) _ ^^(2,2) 


K^^feDgix) = Dg{x) = I {x,y)g'{y)dy 

J a 


d2,l) 

■Af,4 

— J<^(2,1) 


d2,l)/ 


= K)^A{.^,y)g{y)\^- j g{y)—Ky/{x,y)dy. 


Note that 


2Af-l 


d^NAi^^y)- = ^'>Pj{x)^ijkipk{y) 

j,k=0 

2N-1 

= Y1 '^j(^)9jk7^k{y)Vw{y) 


j,k=0 


and 


hence 


K^^’l\x,a) = K)^’^\x,b) = 0 


w{a) = w{h) = 0, 


( 2 , 1 ), 


Continuing, 


/y(2,2) 


eDg{x) 


-/ 9{y)-^K^NA{x,y)dy = J {x,y)g{y)dy 
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i.e., 


N 4 • 


The proof is complete. 


□ 


According to Theorem 13.11 can be written as 


where 


_ 


DK^^fe DK^^’f 




K 


( 2 , 2 ) 

NA 


= : AB, 


A = 



B = 


From fl3.2p and using Theorem 12.61 we have 



/v(2>2) 


K 


K 


( 2 , 2 )' 

V,4 

( 2 , 2 ) 

V,4 


^(2,2) ^(2,2) 

.0-7V,4 ^ ^N,4 

^(2,2) ^(2,2) 


(/)1' = det f / + f as) / ) = det ( / + A (s/) ) = det ( / + (s/) A 


= det J + S(/A 


Since 


B fA) = 


iv(2,2) 

-^'■Af,4 

iv(2,2) 

-^'■Af,4 



^(2,1) 

j^(2,l) 
-^'■Af,4 

-^'■Af,4 ^ 
J^N,4: ^ 


KS:AfD liSff 

Kf^efD 


then 


GfU) 


1 2 


= det 


I + k'-^AI'd 
K'ifefD 


jv(2>2) p 

^N,4 J 

i+ GT/, 


The computation below reduces the above into a determinant of scalar operators. We subtract 
row 1 from row 2, 

Gg)(/)l' = det 


I + K^ifefD K^ff 
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Next, add column 2 to column 1, 




= det 


I + K^^fefD + K^^ff 


^(2,2) r 

-^Ar,4 J 


det (l + 


■fefD + Af;?/) 


Remark. The above result agrees with [6] for the GUE case if we take / = —/i xj, where xj is the 
characteristic function of the interval J. 

Now we use the commutator [D, f] := Df — fD to obtain a better suited result for our purpose. 
For a given function smooth g{x), we have 


[D, f]g{x) = Dfg{x) - fDg{x) 

= {f{x)g{x)y - f{x)g'{x) 
= f'{x)g{x), 


this is, 


lA/! = /'. 


It follows that 


fD = Df-[DJ] = Df-f. 


Taking this into account, we have the following theorem. 


(3.3) 


Theorem 3.2. 


G‘^'(/)] " = det (/ + 2Agf>/ - Agf e/') 


where the kernel of reads 


2N-1 


j,k=0 


3.2 GSE 


In the case of the Gaussian weight w{x) = e , a: G M, we again follow the discussions PE], and 
choose a special ipj to simplify as much as possible. To proceed, let 


1 1 

'<p 2 j+i{x) := -^ip 2 j+i{x), 'Ip 2 j{x) := -^eip 2 j+i{x), j = 0,1, 2,... 


(3.4) 
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where (y2j(x) is given by 


iPj{x) = 


Hj{x) ££ 


", cj = IT* 2" yr(j +1), 


( 3 . 5 ) 


and Hj{x),j = 0, 1 ,... are the nsnal Hermite polynomials with the orthogonality condition 

/ CO 

HHnix^Q dx 

-oo 

^2 

We show in the next lemma that, this dehnition satishes fl3.1l) . i.e., tpjix) = 7 rj{x)e~^, j = 0,1, 2,.. 
where tt^- (a:) is a polynomial of degree j. 

Lemma 3.3. 'ipj{x)e'^, j = 0,1, 2, ... is a polynomial of degree j. 

Proof. If the index is odd, then it is clear that 'ifj2j+i{x)e 2 is a polynomial of degree 2 j + 1 . For 
even index, 

x2 ^2 

^ 2 j(x)e^ = -^{e(p 2 j+i{x))e'^ 


1 

=-^e 2 e 

\/2 

1 


C2j+1 


-e 2 


H2j+i{x)e "2 


H2j+i{y)e ""^dy- / H2j+i{y)e dy 


2\/2c2j-)-i 

-P 

Since H2j+i{y)e 2 ^ j = Q, 1 ,... is an odd fnnction, we have 


J=0,l,2, 


0 = 


H2j+i{y)e '‘■^dy= / H2j+i{y)e dy + / H2j+i{y)e dy. 


Hence 


and we hnd. 


H2j+i{y)e "2 dy = 


f)2j{x)e"^ = 


1 

e 2 


■\/2C2j+l 

From mathematical indnction, it follows that 


H2j+i{y)e “2 dy, 


H2j+i{y)e ''2 dy. 


k-5 


/ y^e ^'2 dy = —e "2 ^— I'jx^ ^ + {k — l){k — 3 )x 

J —OO 

+ ■■■ + {k-l){k- 3 )--- 2 ], A: = 1 , 3 , 5 ,-- - . 

y2 

Since H2j+i{y) is a linear combination of y,y^,..., we see that Id2j+i(y)G~^dy is eqnal 

x2 ^2 

to e~^ mnltiplying a polynomial of degree 2 j. It follows that 'ip2j{x)e~^ is a polynomial of degree 
2 j. The proof is complete. □ 
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Using fl3.4p to compute 


: = 


we obtain the following lemma. 


- 'il)'j{x)i)k{x)) dx 


2N-1 


j,k=0 


Lemma 3.4. (Dieng, Tracy — Widom) 

/ 


= 


0 10 0 

-10 0 0 

0 0 0 1 
0 0-10 


0 o\ 
0 0 
0 0 
0 0 


0 0 0 0 ■■■ 01 

\0 0 0 O--- -10/ 

' / 2Nx2N 

It is clear that ^ = —M^^\ so ii 2 j, 2 j+i = fJ- 2 j+i, 2 j = 1, and /ijfc = 0 for other 

Hence 

2N-1 


j,k=0 

N-l 

= ^2i(x)^/^2j+l(y) - V’2j+l(2^)V’2j(l/) 


N-l 


j=0 

1 
2 


j=0 


'N-l N-l 

‘^2j+i{x)ip2j+iiy) - e<^2j+i{x)Y+iiy) 

.j=0 j=0 


Recall that the Hermite polynomials Hj satisfy the differentiation formulas [10] 

Hj{x) = 2xHj{x) - Hj+i{x), j = 0,1, 2,..., 

H'j{x) = 2jHj_i{x), j = 0,1, 2,.... 

^2 

Using the fact that Hj{x) = Cj ^Pj{x) e^, fl3.7p becomes 




'Ax) = Xipj{x) - \/2(i + l)^j+i{x), j = 0,1, 2,... 


and similarly, fl3.8p becomes, 


A 


Ax) = -Xipj{x) + ^/^ipj-i{x), j = 0,1,2, 


cases. 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 
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(3.11) 


Combining fl3.9p and fl3.10p . to eliminate we obtain 


i=o 

Af-l 




= y - Y 2 + j = 0,1 , 2 ,... 

Using (13.lip to replace ^' 2 j+i{y)) we find, 

Af-l 

^£(p2i+l(x)(pY + l(2/) 

y^J + \y>2j{y) - 3 + l^2j+2{y) 

N-l j - - N-1 

Y-^' + ^£y^2j+i{x)(p2j{y) - + ^£‘f2j+i{x)(p2j+2{y) 

^ r .—1 ^ . I —r 

~ Y1 Vj^V’2j-i{x)^2j{y) - Jn+ -e^2N+i{x)>^2N{y) 

■ - ^ j=0 ^ 

y^e>f2j-ii,x) - \Jj + ]^e^2j+i{x) 


3=0 

N-l 


3=0 

N 

-E 

3=0 


V 23 {.y) - XjN + -eip 2 N+l{,x)ip 2 N{,y)- (3.12) 


To proceed further, using Lemma 12.71 together with (13.lip we find 

ip 2 j{,x) = eD^p 2 j{,x) = e i^' 2 j{x) 

= V~j ^^2j-l{x) - \j3 + \ £y>2j+l{x). 

Substituting fl3.13l) into fl3.12l) . it follows that 

iV-l N j - - 

Y ^^2j+lix)(p2j+liy) = -5^<^2i(a))(p2i(2/) - yA^T -£(^27V+l(x)<^2iv(y). 


(3.13) 


3=0 

Hence, fl3.6p becomes, 

K^N^i^^y) = I 


3=0 


'N-l N I - - 

Y ^ 2 j+iix)^ 2 j+i{y) + Y 7^2i(a^)v?2i(|/) + Y + -e(p 2 N+i{x)(p 2 Niy) 

.3=0 j=0 ’ 


.3=0 

2N 


Y + Y ^ + 2^7^2JV+i(2:)7’27v(i/) 

.3=0 ^ 


SNix, y) + \ N + -eip 2 N+iix)ip 2 N{y) 


where 


Q / \ \ r \ I AT I 1 ‘f2N+l{x)(.p2N{y) - V2N+l{y)‘^2N{x) 

SN{x,y) := 2_^(pj{x)^j{y) = I N + - ---, 

j=o V ^ X y 

and here the last equality comes from the Christoffel-Darboux formula. 

By Theorem 13.21 we have the following theorem. 
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Theorem 3.5. 


^7V^(/) — det / + S'at/ — -S'atS/' + \ I N + - {eip2N+l ® <f2Nf) + o\/^ + O (^V^2iv+1 ® SLP2n) f 


3.3 Large N behavior of the GSE moment generating function 


To proceed with the large N investigation, write, 


C(/) 


, as 


gW(/) =:det(/ + r), 


where 


T Sn/ — -Sn^I' Y 2 ® ^2Nf) 2 Y 2 ® £<d2v) f- (3-14) 


We hnd. 


log det(/ + T) = Tr log(/ + T) = Tr T - ^Tr + ^Tr -. 


This is obtained by the trace-log expansion of logdet(/ -|- u T), where 0 < z/ < 1, follow by a 
continuation to p = 1. We assume that similar continuation also holds in other cases. 


Theorem 3.6. 

1 ^ / X y \ sin(x-y) 

V™ ViiV VViiv’ ~ 7r(x - y) ’ 

,ta 1 

7V^<X>^4/V VVdiV V4iV/ 

The next theorem characterizes the large N asymptotic of various “scaled” quantities. 

Theorem 3.7. 


‘f 2 N 


f—) 

VV4iV/ 


= 


(—l)'^7r 2jV 4 cosx-|-O^A^ 4 j ^ jY ^ oo^ 


£^ 2 N 


~ ^(—l)^7r 4sinx-|-O^A^ 4 j ^ ]Y 00 ^ 




f —1 

VV4iV/ 


-2-iAr-i+ 0 fiV-t ) , N^oo. 
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Proof. From the asymptotic expansion of Hermite polynomials [13], 


HJx) e 2 = 


cos ( \/2n + la:- —^+0[n 2 j ^ qq, 


(3.15) 


where 


r(2?7, + i) r(2?7, + 3) i 

^2n — “iTT— r~rr ^2n+i — ~frr~~n^ + 3) ^ . 


r (n + 1) 


r(n + 2 ) 


We have 


P 2 N 


^_ 

J 7rh^^r(2iV + 1) V\/4iVy 

1 r(2A^ + i) 


(if) 


cos 


7ri2^^r(2iV +1) r(Ar + i) 

1 ^/T{2N + 1) 

714,2,^ r(A^ +1) 

= {— 1 )^ 7 i~^N~^cosx + o(^N~^'^, N ^ 00 , 


AN 

cos(x — Ntt) + O ( N~^ 


^m^,-NA+o{N-i) 


where we have nsed the Stirling’s formnla fl2.ip . 
A straightforward compntation gives, 


e(p 2 N{x) = - 

1 

2 


= 2 '^ 


nx poo \ 

/ (f2N{y)dy- / <^2N{y)dy\ 

' —00 j X / 

^0 px /•O poo 

/ ip2N{,y)dy+ / if2N{y)dy- / ip2N{y)dy - / ip2N{y)dy 
' —00 Jo J X Jo 

nx nO poo 


' P2N{y)dy+ / P2N{y)dy- / y>2N{y)dy 

0 J —00 J 0 

\ f 2 rO 2 r°° 2 

1^,1 W / H2N{y)e~^dy+ H2N{y)e~'^dy- H2N{y)e~^dy 

7r42^+Vr(2iV + 1) V Jo J -00 Jo 

I 2 

H 2 N{y)e~^dy. 


7r42^^r(2Ar + 1) Jo 

So we hnd, for large A^, 

(*) 


VIn 


H 2 N{y)e "2 dy 


7r42^^r(2iV + 1) io 
—i- , = • , / i?2V ( —:== 1 e~^dy 

7r42^^r(2iV + 1) \/4]V Jo V\/4iV/ 

1 1 r(2A^ + i) r 

7rh^^T{2N + 1) y/AN r(Ar + l) Jo 
y/n 2 N + l) 


(—1)^ cos y + O [ N 2 ) dy 


1 ^_ (—1) sinx + O 2 

7r4 2^^/4iV^(A^ +1) L' V 

= 2“^(—l)^7r“2iV“4 sinX + O ^A^“4 j ^ ]Y 00 ^ 
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where we have used the Stirling’s formula fl2.ip in the last step. 
Finally, a straightforward computation gives, 

" rx roct 

e<^2N+i{x) = 9 / ^2N+i{y)dy- / (p2N+i{y)dy 


_J —oo 


\ 2 2 
1 ^,3 / H 2 N+i{y)e~'^dy - / H 2 N+i{y)e~'^dy 

7ii2^+2^T{2N + 2) U-oc Jx 

^ r 2 2 

/ Ff2V+i(2/)e“^d2/+ / //2Ar+i(y)e-^d2/ 

7r4 2 ^+ 2 ^r( 2 A^ + 2) U-oo io 


) 2 2 
H2N+i{y)e~^dy - / i/27v+i(?/)e“^d?/ 


i 

7r32^+i^r(2A^ + 2) 




2 2 
^2V+i(2/)e“^d|/- 2 / hf27v+i(|/)e“^dy 


continuing, we see that. 


r 2 /■°o 2 

.i2"nyr(2jv + 2) U -y„ 


£y^2v+i 


a: \ 

.71^; 


By fl3.15|) . we have 


.i2-H jr(2N+2) Ir 

1 f°° _jd 

.i2«nvr(2^ + 2) io 

7rh«+4 yr(W+l' I vi^)" “ 

/•°0 2 


FT. 


2Ar+l 






OO. 


H. 


It follows that 

(tI^) ^ ^tIn + 2 ) ^ (-1)^(1 - cosx) + O [n-^) 

On the other hand, from [S], 

. s _jd 2 (-l)^r(2A^ + 2) / 3 \ 

I -fyvTT)-2np-vi;5;2j. 

It follows from Lemma [2.21 that, 

I H2N^dy)e-‘^dv=^‘- 


, N ^ OO. 
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By the Stirling’s formula fl2.ip . we obtain 


^H^2N+1 


= 2-^(-l)^7r-5Ar-i(i-cosx)+ 0 - 2-5A^-i + O 

= —2~^N~i+o(^N~i^, N ^ oo. 


□ 

We are now in a position to compute TrT and TrT^ as A^ —)■ oo, using Theorem [XS] and Theorem 
13.71 The estimates provided by Theorem 13.71 are instrumental in the large N computations that 
follows. 

In what follows, we replace f{x) by / ( y/ANx ), and note that 


The f that appears in the trace will be accordingly interpreted. 
We hrst consider Tr T, which reads, 


Tr T = Tr S'at/ - Tr ^SnsJ' + Tr ^A^ + ^ {e<f 2 N+i ® 7 ^ 2 ^/) + Tr N + ^{eip 2 N+i ® £^ 2 N)f- 


Tr T has four parts. 
First of all, we hnd 


TrS'jv/ = f S]s[{x,x)f (^VaNx^ dx 


' —OO 
bOO 


-.s 


J-oo ^/4iV"^ V\/4iv’ ^/4]V 

1 

—)■ — / f{x)dx, N ^ oo. 


f- 


X x 


f{x)dx 


The second term reads. 


= 1 


poo poo 


SN{x,y)e{y,x)f ('/ANx) dxdy 


1 7“ 


— OO J —OO 
POO 


—oo L»/ X 


SN{x,y)dy 


SN{x,y)dy 


f ^\/4A^xj dx. 


To proceed further, let 


u 


= VANx, V = y/ANy, 
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it follows that, 


1 1 r 

Tr-SNef = 

2 ^ AV^J-c 


-.S 


V\/4iV’ y/ANj 

■dv — i 


(- 


u V 


-)■ 


1 


4\/4iV7r j-o, 

Now let t = M — it follows that, 

Tr^SNsf ^ 


I dv — 


-.S 


u V 


-oo V\/ 4 iv’ y/m) 


dv 


f{u)du 


sin(M — v) sin(tt — v) 


dv 


u — V 


u — V 


f{u)du, N ^ oo. 


Ay/ANn J- 
1 


sinf , smt ' 

- dt - / - dt 


OO —oo 
oo 


f\u)du 


A^/W^lJ-oo ^2 

0, N ^ oo. 


^ - I) f'{u)du 


For the third term, we have. 


Tr )Jn {eip 2 N+i < 8 ) <f 2 Nf) = ^ J £^2N+i{x)(p2N{x)f (yiNx^ dx 


= \lN+l-—^f eip2N+l (^2N ( ^ 


2 y/m J- 

/-oo 


2V27iN J- 

where use has been made of Theorem 13.71 

Finally to the fourth term, and take note of Theorem 13.71 


cosx f{x)dx + O , N ^ oo, 


f{x)dx 


N +-{e(p2N+i® £V2N)f — 2\ ^ ^ 2 J 


- Viv+i.' r 

2V 2 y/mj-oc 

= O (N-^) , N ^oo. 


£V^2N+i{x)eip2N{x)f \^\/ANxj dx 

( X 


f;®) (;lf) 


Therefore, the large N expansion of TrT, reads. 


TrT = 


TT 


f{x)dx — 


- 1 ) 


N roo 


2y/27iN J-c 


cosX f{x)dx + O (^N , N ^ oo. 
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Working out Tr with T given by, fl3.14p . there are 10 traces: 


TrT^ = Tr S ^fS^f - Tr S^fS^^ef + Tr y/WT 2 SNf{e^ 2 N+i ^ ^2Nf) 

+ Tr \jN + ]^SNf{eip2N+i ® £^2N)f + Tr ^SNsfSNef 

- Tr )JN + ^SN£f{eip2N+i O 9^2^/) - Tr N + ^SNef{e(p2N+i <8 £9^2^)/' 


+ Tr + -j {e(p 2 N+i O 9^27v/)(£952v+i O >f 2 Nf) 

+ Tr + 0 (£952Af+l ® 9 ^ 2 Ar/)(£<^ 2 Ar+l ® eip 2 N)f' 

+ I (£9^2Af+l ® £^2N)f'{£^2N+l ® eip2N)f- 


In the following, we calculate the trace on the right side of fl3.16p . term by term. 
The hrst term: 


TiSn/SnI = 


roo roo 


’ —00 ^ —00 
poo poo 


' —00 J —00 


SN{x,y)f (yiNy^ SN{y,x)f (yiNx^ dydx 


-.S 


V\/4iV’ y/mj 


( 


X y 




f{x)f{y)dxdy 


POO POO 

sin(x — y) 

7 r 2 / / 

J —00 J —00 

1 - 

1 

1_ 


f{x)f{y)dxdy, N ^ 00 . 


The second term reads, 

TiS^fS^ef = 


00 poo poo 


S]\f{x, y)f (yj4:Ny') S^iy, z)e{z, x)f' (\/ANx \ dxdydz 


— 00 J —00 J —00 
poo poo 


( 3 . 16 ) 


Sk(.t, !/) 


' —00 J —00 


SN{y,z)dz 


SN{y,z)dz 


f ^\/4iVa;j / dxdy. 


A change of variables. 


u 


= \/ANx, V = V4:Ny, w = y/ANz, 


give 


TrS^fS^ef = 


r*oo poo 


1 / u 

= Ojv 


(- 


4 


2y/AN J-ooJ -00 \y/AN’y/ANj 

(dw- ^ 


I ^ f V 
-Z>N 


-9 


y/AN \y/m'y/ANj 

1 /■°° sinfa - v) 


-.S 


-00 y/ 4 N^^ V\/ 4 iV’ y/ANj 


( 


V w 


] 


dw 


f'{u)f{v)dudv 


2y/4N7i^ J -00 J-oc u-v 
f'{u)f{v)dudv, iV —)■ cx). 


sin(n — w) 


V — w 


dw 


sm n — w] 


V — w 


dw 
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To proceed further, let t = v — w, we see that, 


sin u — w 


dw — 


sin u — w] 


dw = 


-dt = 2 


= 2 Si(n —-u), 


where Si(a:) is the sine integral 


Since Si(—x) = —Si(x), it follows that 


Si(^) r 


TtS^fSp/ef ^ 


Sin n — V 


J-ocJ-oc u-v 


sin -u — V) 


J-c 


Si{v — u)f{u)f{v)dudv 
— Si(n — v) f (u)f {v)dudv, iV —)■ cx). 


The third term: 


Tr y/AN + 2S'7v/(£V?2Ar+l <8 V?27v/) 


= VAN+ 2 


SN{x,y)f (y/ANy) e(f2N+i{y)‘A^2N{x)f (y/Wx) dxdy 


= (tif’tif) [Am) 


y/AN 


(-ir r 

7rV27iN J- 


N poo poo 


sin(x — y) 
x-y 


cosx f{x)f{y)dxdy + O [N ^) , N ^ oo. 


The fourth term: 


Tr W + -SNf{e(p2N+i 0 eip2N)f 


= 


SN{x,y)f (y/ANy) e(p2N+i{y)£(P2N{x)f (y/^x) dxdy 


2 ' Xw i-ooi-oo ( VIw’ J^n) ( viw) ( VIw) 

= O (N-^) , N ^oo. 


X y 


The hfth term, with the change of variables. 


= VANx, V = VANy, w = VANz, r = VANt, 
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we see that, 


Tr ^S^efS^ef 


^oo poo 


64N 


' —OO J —OO 


1 / u 

iOtv 


I I dv — 


1 c- ^ ^ 

-ON 


(- 


= O (N-^) , N ^oo. 


y/4N \y/4N'y/4Nj 

“ 1 C, ^ ^ 

ON 


1 / u 

zOn 


(. 


" Idr- 


( 


-OO V4iv VV4iV’V4iV 


dv 


" ^dr 


-ooViiV \y/4N’y/4Nj 


f'{u)f'{w)dudw 


The sixth term, 


Tr \j^ + -;^SN£f'{s^ 2 N-^i ® ^ 2 Nf) 


N + - 
8N\ 2 


r*C 30 poo 


—OO J —OO 


^ 2 N 


U \ / W 


VIn) 


£^2N+1 


ViiV 


1 o _ 

ViN \vm' vmj 

f iy) f {w)dudw 


^ (it; — 


—'l 

./m ” VV4/v’ Vt/v/ 


1 c- ^ “ 

-ON 


dv 


= O {N-^) , N ^oo. 
The seventh term, becomes. 


2 V ^‘p2N)f' 


16iV V 2 


£^ 2 N 


(*oo poo 


— OO J —OO 


1 / u 

zOn 


( 


y/AN VVT/V’VI/V/ 


I (in — 


1 / u 

zOn 


( 




-ooVan VVt/v’Vt/v/ 


dv 


(A)-™'-' 

= O , N ^ 00 . 


£<-P2N+i J f'{u)f'{w)dudw 


The eighth term can be computed in a similar manner, and we have. 


Tr ( iV + - ) {eif 2 N+l ® ^ 2 Nf){.eV 2 N+l ® ^ 2 Nf) 


foo M \ 


AN 

= O (N-^) , N ^ 00 . 


(^ 2 N (^ 2 N ( - 


£V2N+1 I r— ) £V^2Ar+l I I ‘f2N I r— V2N 1 r— 

-ooJ-00 \\/AN J \y/AN J \\/AN J \\/AN 


f{u)f{v)dudv 


Proceeding in a similar manner with the ninth term, we have. 


Tr (£V? 2 Ar+l ® V? 2 Ar/)(£<P 2 Ar+l <8) e<P 2 Ar)/' 


N /-OO .00 


(- 


u 


4iv y-ooloo^^'^^^viiv 


£'^ 2 N +1 


Vm 


) (tw) (tto) f>Wv)d'udv 


= O , N ^ 00 . 
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The tenth and last term in TrT^, becomes 

Tr ^ + 0 {e^p 2 N+l ® £^ 2 N)f\s‘^ 2 N+l ® £^ 2 N)f 


N + ^ f°^ f°° 


—oo ^ —oo 


(;^) (;w) (;w) (;w) 


16A^ 

= 0{N-^) , N ^oo. 

Hence, the large N behavior of fl3.16p . reads, 


TrT^ = 


TT^ 


roo poo 

sin(a: — y) 

J —oo J —oo 

x-y 


f*00 /‘OO 


sin(a: — y) 
27r2x/]V 7-00 7-00 x-y 
(—1)^ f°° f°°sin{x — y) 


f{x)f{y)dxdy 
Si(a; - y)f\x)f\y)dxdy 
cosx f {x) f {y)dxdy + O {N~^') , N ^ oo. 


7 rV2TiN 7-00 7-00 x-y 
We are now in a position to compnte the mean and variance of the (scaled) linear statistics 
F ( VANxj ), which are obtained as the coefficients of A and A^, of logdet(J + T). 


Since 


/ (VAN: 


X 


\ 2 

-XF (VINx^ + yF^ (^\/4iVx) , 


we replace / with —XF + in the expression of Tr T and Tr T'^. A minor rearrangement gives. 


log det (/ + T) 


N /-oo 


= -X{ - r F{x)dx - ^ / 

vr7_ 2\/2tiN J-o 


— OO 

2 11 Foo 


A 

^ 1 TT 
1 


cosx F(x)dx + O (A^ 


F{x)F{y)dxdy 


— OO 

poo poo 


2 poo roo 

sin(x — y) 

7r2 / / 

J —OO J —OO 

x-y 


271^VN 

^ ^ / cos X F^(x)dx - 

2V2^ [J-oc tt 


OO j —OO ^ y 

N r 


sin(a; y) ^ y)F'{x)F{y)dxdy 


f*oo poo 


sin(x — y) 


cos X F{x)F{y)dxdy 


+ 0 {N-^) I, N^oo. 


I—OO J —OO ^ y 

Denote by and the mean and variance of the linear statistics ^ ^Nxj^, then 

we have obtained, the large N corrections of these qnantities. 


Theorem 3.8. As N ^ oo, 


(GSE) 1 (GUE) 

hiv = 2 ^^ 


-1) 


N poo 


4\/27rA^ J-c 


cos X F{x)dx + O [N ^) , 
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V 


(GSE) 

N 


= -V 


(GUE) 

N 


f*00 POO 


2 ■ 47 ^ 2 ^/]V J- 


Si(x - y)F'{x)F{y)dxdy 


OO J —OO 


X - 2/ 


TV 


^ ^ / COS X F^(x)dx - 

4\/M Li-oo tt 


r*C50 /*oo 


sin(x — y) 


cosx F{x)F{y)dxdy 


' —OO J—OO ^ y 

where and /or N ^ oo are given in / li.^l) and U.3\} respectively. 


+ 0{N-'), 


3.4 LSE 

We study the case with the Laguerre background, namely, the weight, 

w{x) = x"e“*, a > 0, x G M+. The idea is to choose special so that takes on the simplest 
possible form. To this end, let 

1 


'ip2j+i{x) := ^(pg+i^(x), / = 0,1, 2,... 


i^2jix) := -^£V?^“+i^^(x), j = 0,1,2,... 


where (x) and (x) are given by 


‘h’j (^) = (a-1) ^ J =0,1,2,..., 

r • 

^3 


•~(q;— 1) / \ 

w = 


S._l _£ 

d O J- y-v O 




I! e ", i = 0,1,2,... 


(3.17) 

(3.18) 


(3.19) 


Here L, (x). j — 0,1,2,..., are the Laguerre polynomials, with the orthogonality condition, 

^ L^^\x)Ll3\x)x^e-^dx = Sjk, cf = 


It is easy to see that 


^^(x)(pj," ^\x)dx = 6jk, j,k = 0,1,2, 


We now prove that (13.171) and fl3.18p satisfy (13.ip . i.e., i^j{x) = Ej{x)x 2 e 2 , j = 0,1, 2,..., where 
Ej{x) is a polynomial of degree j. 


Theorem 3.9. 'ipj{x) = Ej{x)x‘^c 2 , j = 0,1, 2,..., where Ej{x) is a polynomial of degree j. 
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Proof. We prove this by considering two cases, j odd, and j even. If j = 2?7, + 1, then by fl3.17p . 




\/2c. 


x). 


2n+l 


Let j = 2n, then 


'lf2n{x) = - ^e(p2n+l{x) = - ! [ 4n+?(2/)l/^ 

\/2 V 24 „+/ io 

and we have nsed the fact {y)y^~^e~^dy = 0, n = 0,1, 2,.... 

Let ns rewrite the above as 


L'fn+\\y)y'^ ^^dy = iT 2 n{x)x'^e 2 , n = 0,l,2,... 


a X 


(3.20) 


where 


Take a derivative on both sides. 


^2nix) := -V2C2n+^^TT2n(x). 


L’^+^i{x)x2 ^e ^ = (^XT2n(x) + ^^2n{x) - ^X92n(x)^ X2 ^e 2, 


which becomes. 


Lm{x) = X7r'.^^{x) + ^^2n{x) - ]^Xn2n{x). 


Note that equation fl3.2Up is eqnivalent to fl3.2ip . Now we seek to solve fl3.2ip . 
Snppose 


(3.21) 


^2n{,x) := ao{n) + ai(n)x + a2{n)x H-h a2n-i{n)x + a2n{n)x , 


we see that the right side of fl3.2ip is eqnal to 


tt 1 

X%nix) + ■^^2n(x) - -XW2n(x) 


fao(n) 


(l + I) ai(n) - ^ao(n) 


X + 


-f 


n 


ax 1 

2n - 1 + - j a 2 n-i(n) - -a 2 n- 2 (n) 
^a2n(^^)a;^’"■^^ 


+ 


a 

2?7, H— 
2 


-a\{n) x^ + 


n) - -a 2 n-i{n) 


X 


2n 


(3.22) 
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On the other hand, the left side of fl3.2ip is equal to 


r(“-l) 

-^2n+l 


(x) 


f2n + a\ 

V2«+V 

1 

'"^W) 


(2n + a\ I (2n + a\^ _^(2n + a\ 

V 2n ) 2\\2n-l) ?,\\2n-2) 

2n + a\ 2 n _ 1 (2n + a\ o^+i 

I ) (2n + l)!V 0 ; 


X 


Compare the coefficients of fl3.22p and fl3.23p . we have the equations 


(l + f)o,(n)-P„W = -(“r) 

(2 + f) 02(n) - pi(n) = 

(3 + I) a^(n) - ^a^{n) ^ 

< 

(2n - 1 + f) a 2 n_i(n) - ^a 2 n- 2 in) = 
(2n + f) a 2 n{n) - ^a 2 n-i{n) = 

1 / \ 1 /2nH-Q\ 

~ 2 ^‘in{n) — — (2n+l)\ ( 0 )' 

By solving the first 2n + 1 equations in fl3.24p . we find, 


(3.23) 


(3.24) 


i“o(«) = hS?) 

h.(")=- diOr) 

1 ( \ ^ 1 /2n+a\ _ 1 (2n+a\ , 1 1 (2n+a\ 

a{a+2){a+i)\2n+l) (o+2)(a+4) V 2n / 2! o+4 V2n-1/ 

< (ri'\ _ _1_ /2n+a\ _ 1 _ /2n+a\ . 1 1 /2n+a\ _ 1 1 /2n+a\ 

2^3\iI'} — a{a+2)(a+4:){a+6)\2n+l) (o+2)(Q;+4)(a+6) V 2n / 2! (o+4)(a+6) V2n-1/ 3! a+6 V2n-2/ 


1 (rt'] = _1_ (2n+a\ _1_/2n+a\ , J__1_ /2n+a\ 

2^2n\ ) 0(0+2)(a+4)"-(Q:+4n) \2n+l/ (o+2)(o+4)---(o+4n) V 2n ) ' 2! (o+4)(o+6)"-(o+4n) \2n—1/ 

_J__1_|'2n+o\ I . . _1_1_/2n+o\ , 1 1 /2n+o\ 

3! (o+6)---(o+4n) V2n—2/ (2n—1)! (o+4n—2)(a+4n,) V 2 / (2n)! a+4n V 1 /' 

The last equation of fl3.24p . simplifies to. 


a2n(,n) 


2 

(2n + 1)!‘ 


By Lemma [2.31 we see that linear system fl3.24p is solvable. Hence '^ 2 n{x) is a polynomial of degree 
2n. The proof is complete. □ 





































With fl3.17p and fid.isp . we compute := — 7/jj(x)ipk(x)) dx).^^^, resulting 

in the following theorem. 


Theorem 3.10. 


= 


^0100 
-10 0 0 

0 0 0 1 

0 0-10 


0 0^ 
0 0 
0 0 
0 0 


2Nx2N 


0 0 0 0 ■■■ 01 

\o 0 0 o--- -loy 

Proof. Let rrijk be the (j, fc)-entry of i.e., 

poo 

rrijk := / {'iljj{x)'i/j',.{x) -'i/jj{x)'ipk{x)) dx, j, k = 0,1,... ,2N - 1. 

Jo 

We compute by considering four cases: (j,/c) = (even, odd), (odd, even), (even, even), (odd, odd). 
For the (even, odd) case, 

poo 

m2j,2k+l = / {'i/j2jix)iJ2k+i{x) - ir'.2ji^)'^2k+i{x)) dx 

Jo 

poo poo 

= / ij2j{x)'ljj2k+lix)dx - ij2j{x)fj2k+l{x)dx 

Jo Jo 

( poo \ poo 

fr 2 j{x)ij 2 k+l{x)\^ - J i> 2 j{x)'llj 2 k+l{x)dx] - J ij 2 j{x)'ll 22 k+l{x)dx 

poo 

= -2 J 'il;2j{x)'ilj2k+i{x)dx 


= -2 






(pij+l\x)cp^2k+l\^)d^ 


= d. 


jk- 


For the (odd, even) case, since is antisymmetric, we have 

m2j+l,2k = —m2k,2j+l 

— ~dj^k- 
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For the (even, even) case and j > k, 


m2j,2k = / {'Ip2j{x)ij2k{x) - 'll^2j{x)lp2k{x)) dx 

Jo 

poo poo 

= / 'Ip2j{x)i>2k{x)dx - / 'ijjL{x)'tp2k{x)dx 

Jo Jo 

/ OO 

1p2j{x)ll)2k{x)dx 

poo 

= -2 / 'ilj'{x)'ip2k{x)dx 


= -2 


—-j^L^^j+l\x)x"^ 2 j ( 7 r 2 fc(a;)x 2 e 2) dx 

V 2C2,-I 1 / 


/ L^^j+l\x)n2k{x)x'^ ^e ^dx 


^2j+l “^0 


= 0 , 


since 'K2k{x) is a polynomial of degree 2 k which is less than 2 j + 1 . 

For the (odd, odd) case and j > k, 

poo 

m2j+l,2k+l = / {'^|J2j+l{x)^p2k+l{x) - 'llj2j+i{x)'ljj2k+l{x)) dx 

Jo 

poo poo 

= / J^2j+l{x)ij2k+l{x)dx - / ^l)2j+l{x)'^jJ2k+l{x)dx 

Jo Jo 

/ OO / poo 

J^2j+l{x)J^2k+l{x)dx-i'ljj2j+l{x)'ljj2k+l{x)\’^- / 'ljj2j+l{x)'llj2k+i{x)dx 


= 2 / ij2j+i{x)ij2k+i{x)dx 

Jo 

= (^-1)^ (^-1) [ [x te;;^(a:)) + |4“;;^(a;) - ^a;4“;;^(a;)') x^-V^da; 

C2j+ic^2k+i Jo V ^ ^ ^ 


(^-1) (^-1) / 4“+!V) (a; (4fc+?(a^)) + 14fc+l^(a^) - j x“ ^e "dx 

%'+l '^2fc+l “'0 V / 


= 0 , 


since x ^L2'^_,_J^^(x) j + (x) is a polynomial of degree 2 A: + 2 which is less than 


2j + 1 . 


If j < k, due to the fact that is antisymmetric. 


'm-2j,2k — —Tn2k,2j 


= 0 , 
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^2j+l,2fc+l 


—fn2k+l,2j+l 


= 0 . 


Thus 

m2j,2k = m2j+i,2k+i = 0, j, /c = 0,1,..., - 1. 

This is just the desired form of 

It’s clear that ^ so H 2 j, 2 j+i = — l,h2j+i,2j = 1, and fijk = 0 for other 

The rest of this subsection is devoted to the determination of y), 


2N-1 




j,k =0 

N-l N-1 

^2i(a^)^2i+i(2/) - Y ^2i+i(a;)^2i(2/) 

j =0 j =0 

^ N-l N-l 

- Y -oY \^%'+hy) 


3=0 


3=0 


From fl3.19jl . we see that, 


(a—1)/ \ 
Lf) ’(X) 


d“-l) 


^ _ -| _ ^ 

x 2 e 2 


Recall that the Laguerre polynomials ^^(3^)) satisfy the differentiation formulas |8], 

X (^Lf~^\x)^ = iLf~^\x) -{i + a- l)Lfs^^\x), j = 0,1, 2,..., 

X = U + ^)lYi^\x) + {x-j- a)LY^\x), j = 0,1,2,.... 

Summing fl3.26p and fl3.27p . and divide by 2, gives, 

X fL<“-‘'(x)V = J = 0.1.2.,, 


x{lY"\x))' + Yt^lY'\x) 


or 


X I Li"-‘'(x) I + ^L';-‘'(x) = ^L-7"(x 


a — X 


(«-!)/ 


i + 1 r(»-l)^ 'i j + a — 1 ^ (g-l) 


L^^_r\x), j =0,1,2,.. 


Hence fl3.25p becomes. 


(a-l)/ \ 

ip) ’[x) 


d“-l) 


^_1 _ X 

X 2 e 2 


LtPix) - l±^Xl&-p{x) 


2 




□ 

cases. 


(3.25) 

(3.26) 

(3.27) 


(3.28) 


31 

































We see that 


^3 


(q-1) 


is a linear combination of and just like the GSE case studied 


(a-l) 


in the last section. Replacing j by 2j + 1 in (I3.28p . to hnd, 


{y) = \l U + 1 ) (j 1 ^2^+2 ^) U + i) ‘^'2“ i.y) 


a; + 1\ _(q_i) 


^ \ ~(a—l) j 


A straightforward computation shows that, 

N-l 


H^^2j+^hy) 


j=0 


X] (^) \J(y + 1 ) (3 + ~ ^ (^)\j + 0 {3 + f) (2/) 





From fl3.28p and Theorem 12.71 


la—11/ \ 

W) 


= e 


= £ 




^ V(2i + l)(2i + a)cp^"+iHa;) - ^ V2j(2i + « - 

i + (i + - \Ij (i + 


Hence, the sum, (^) 7’2?+i'^(2/) 




simplihes immediately, and leads to. 


N-l 


N 




i=o 




-(«-!)/ 


i=o 


a' 

2, 




( 2 /)- 
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It follows that, 


JV-l N 

y) = 7 ; Y1 + oJ2 


j=0 


j=0 


= ^<5jv(a;, 2 /) - f) ^^27+i^3;)^^7^^(2/), 


where 


2Af 


^iv(x,|/) := ^^(a:)v?S“ ^^( 2 /) = -^/{2N + 1){2N + a) 

j=0 


^^^N+li^)A ^\y) - ‘fm+iiy)y^m ^’i^) 




x-y 


Here we used the Christoffel-Darboux formula in the last equality. 
By Theorem 13.21 we have the following theorem. 



3.5 Large N behavior of the LSE moment generating function 


Now consider the scaling limit of 



write 


afif) 


det(/ + T), 


where 
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( 2 2 \ 

8N'8NJ ^ 


where B^°‘ ^\x,y) and B^°‘ ^)(x,x) are given by M.Gi) and \1. ?| j with a replaced by a 

Theorem 3.13. 

y O \ T / \ 

+ 0(^A^-i), N^oo, 


v;r« " ( 1 ^) = 2 ( 2 A')> 


‘J'/’lw '* ( gjv) ^ (2^^) ’ 


Ja-i{y)dy - 1 


L~'o 


+ O (^N 2 j , N ^ oo, 


ef-'S'tl (^) = <2'^)^’ [ + O (jv-i) , N^x,. 


Proof. Recall the asymptotic formula of the Laguerre polynomials 

■ J„ (y(4JV + 2a + 2)x)+xiO (iVt-5) , 


L'P{x)x^e-i = fiV + “ + 


r(iv + i) 


We hnd, 


-(a-i) I a; 


'r( 2 iv + i) 


X 


r-1 


r(2A^ + a) V8^y V8^ 

'r(2iv + i) / x^ rr(2iv + a 


e i6iv 


r(2A^ + a) \8N 

,1 Ja—l(^X^ 


r(2iv +1) + i) ' -^o-iW + o 


= 2 {2N)^ +0 , A^^oo, 


where we have used the formula, 


r(^ + a) _ ^a-b 
T{n + b) 


[l + O (u ^)] , u -)■ CXD. 


Proceeding to we have, 


ep\N (3^) 


Pm^\y)dy- I Pm'’iy)dy 


(a-l) I 


UO 


1 /r( 2 iv + i) 
2\/ r(2Ar + a) _ 


1 /r( 2 iv + i) ; 
2y r(2Ar + a) f 

1 /r( 2 iv + 1 ) 
2V r(2Ar + a) 


^oo 

/ L^°~^\y)y^-^e-Uy 

J X 

f*X poo 


2y^ L<r>(,),»-e-t,, 


2fr(iV + |) 
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N ^ oo. 



































where we have used the fact that 


2 fr(iv + f) 

\y)y^ e 2dy=—-^^—-^. 


r(A^ + i) 


Continuing, 


~(o—1) 


1 /r(2A^ + i) 1 

2\ T{2N + a) m 


(a-i) ( y^\ ( y^ 


8N \8N 


^ 2fr(Ar + f) 
^ ^ r(iv + i) 


= {2N)-2(2N + 


ay^ /r(2iV + a) T 
2 J V r(2iV + l) io 


Ja-i{y)dy - 22 


_i /r( 2 iv + i) r(iv + f) 
V r( 2 iv + a) ■ r(iv + i) 


+ 0 A^-5 


= (2iV)-^ ^ J„_i( 2 /)d|/-l +o(iV-t), iV^oo. 


Similarly, 






1 / r(2iV + 2) 

2y r(2A^ + a + l) _ 

1 / r( 2 iv + 2 ) ■ 
2y r(2A^ + a + l) _ 

/ r( 2 iv + 2 ) r 
\j T{2N + a + l) Jo 

where we have used the fact that 


2 c?^- I ^e “^dy 


( 0 - 1 )/ N 

^ ^ Ol\oi O -*-<3 O / 


2 / 2(i|/- / L^“ ^(l/)l/2 ^e 2(iy 




^2Vyi(l/)l/" 2(iy, 


4vyJ(2/)2/^ 2(iy = o. 


It follows that 


~(a-i) / a;- 
^^2iV+l ( 


T{2N + a + l)Jo ^ 27 v+i(I/)I/ e dy 


^ I r( 2 iv + 2 ) [%(»-!) fiT] 

AN\J T{2N + a + l) Jo \8N) \8N) 

= {2N)-^ f J«_i(?/)d|/ + 0(iV-i), At ^ 00 . 


e-^ydy 

a AT \ a AT 


J^.i{y)dy + 0(N-2 
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□ 


We now use Theorem 13.121 and Theorem 13.131 to compute Tr T and Tr as —)■ cxd. In the 

computations below, we replace f{x) by / (^VSNx^ and f'{x) by 

d 




Consider Tr T, which reads, 

TtT = S„/ - TV ^Sr,Ef' - TV d(^Ar + 1 j (at + |) 0 

- Tr 5 ^ + 5) (^ + f) ® e*?™‘’) /'■ 

So we compute Tr T by calculating the four terms in the right side. The hrst term, 

Tr SnI = J S'Ar(x,x)/ (^VSNx^ dx 

= r —Sn f{x)dx 

Jo AN \8N' 8NJ ^ 

poo 

B^'^~^\x,x)f{x)dx, N ^ 00 . 


-)■ 


The second term, 

Tr^S'jvs/' = ^ Sn{x, y)e{y, x)f' 8iVxj dxdy 


1 


Let 


then 


1 1 r 

Tr-SNsf = - 

2 ■' 16iV Jo 


I SN{x,y)dy- / SN{x,y)dy 

0 IJx Jo 


u = VSNx, V = \/8Ny, 


—Sn (—, —^ vdv - [ —Sn ( —, — 1 vdv 
AN \8N'8N J Jo AN \8N'8NJ 


f ( V8Nx) dx. 


,2 ,,2 \ 




16iV 


^\u,v)dv — / ^\u,v)dv 

Jo 


u f{u)du 
u f{u)du, N ^ 00 . 
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The third term, 


tm/® 


'^+i) 


£</^2'n+J(x)<^2‘^ ^\x)f (ySNx^ dx 


(A^+ 1 ) (A^+ “) 1^00 f x^\ _(c,-i) f x'^\ 

- 4 iv- I [-^)^f(^)dx 

poo r pX 

/ / Ja-i{y)dy Ja-i{x)f{x)dx + 0 [N~‘^^, N ^ oo. 


The fourth term, 


\ \/(^^ ^) f) [ ^^ 2 Ar+i(3^)/' (VSA^x) dx 


> + i) (iV + I) 


>-i) y a;2 \ („_i) / x2 


^V^2Ar+l 1 gjY ) ^‘d2Af 


X f\x)dx 


Jc,-i{y)dy / J„_i(2/)d|/- 1 x/'(x)dx + O (iV ^) , N ^ 00. 


Therefore, 


TrT = / 5(“-^yx,x)/(x)dx- 


Ja-i{y)dy J„_i(x)/(x)dx 


16-^ Jo Ux 
1 /■“ r p 


fi*-" ^\x,y)dy — / s’-" ^\x,y)dy x f'{x)dx 


Ja-i{y)dy / Ja_i(?/)d|/- 1 X/'(x)dx + O (iV ^) , N ^ 00. 
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Next, we compute Tr where there are 10 traces, 

TrT^ = Tr S^fSr^f - Tr Sr,fSr,ef - Tr 2^(^iV + 0 (iV + |) S^f 0 

- Tr ^(^+0 {N+^)s^f (£^( 7 ^ 1 ; 0 £^ 7 -')) f + Tr ^S^efS^ef 

+ Tr + 0 (^ + f) ® ^2‘^“^V) 

+ Tr ^ ^(^iV+ 0 (^ + f) ® /' 

+ Tr (^iV + 0 + I) 0 ^ 2 r^V) (^^27+J ® 

+ Tr (^iV + 0 (^N + I) (£^ 2 ^+] 0 ® ^^ 2 ^^^) f 

+ Tr ^ + 0 (^iV + I) (e^^7+] 0 f ® /'• 


In the following, we need to calculate the traces on the right side term by term. The hrst term. 


TrS^fS^f = 


SNix,y)f (^/SNy) SN{y,x)f (VSNx) dxdy 


rr>iis4)a^-{&4)’‘-<''«“' 

/ OO POO 

/ B^'^-^\x,y)B^'^-^\y,x)f{x)f{y)dxdy, N ^ 00 . 


The second term. 


TrSNfSNef' 




4Ar V8A^’ m 


O AT I I A AT (oAr’OAT- / A AT ( O AT O AT I A AT ( O AT O AT 




AN \8N' 8N 




AN \8N' 8N 


u f'{u)f{v)dudv 


^\u,v) / S*'" ^\v,w)dw — / i?*-" ^\v,w)dw u f'{u)f{v)dudv, N ^ 00 


The third term. 
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The fourth term, 



8iV’ 8N 




8n) 


~(o—1) 


m 


u f'{u)dudv 


Ja-l{t)dt 


u f'{u)f{v)dudv + O [N ^) , N ^ oo. 


Tr -S^efS^ef 


S]\f{z,t)e{t, x)f (^VSNx^ dxdydzdt 



SN{z,t)dt — / SN{z,t)dt 

Jo 


vJ \ 


r 1 

/o 4N^^'\8N'8NJ 


rdr 


u w f {u) f {w)dudw 





The sixth term, 




+ O (N-^) , N ^oo. 


{u,v)dv 


^\u,v)dv 


Ja-l{t)dt 


J L^o 


Ja-l{u)w f'{w)f{u)dudw 
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The seventh term, 



The eighth term, 


Tr + - 

' 2 


N+- 


r(«-l) 


(a-1) 


^‘^2N+1 ® 


~(a—1) „ ~(o—1) 
^^2N+1 ® V^2Af 


(iV + i) (iV + l) 


oo /*oo 


16iV2 


OO poo 



-1) ^ \ f \ f \ 


Aa-l) 


(a-1) 


4a-l) 


0 JO 





0 ^0 


Ja-l(t)dt 


\8NJ 
Ja-l{t)dt 


\8NJ 


\8NJ 




8N J 


Ja-l{u)Ja-l{v)f{u)f{v)dudv + O {N ^) , 


The ninth term. 


uv f {u) f {v)dudv 
N ^ oo. 


Tr + 0 (^N + I) 0 ^V) (^<^ 2 '^+] ® f 

(iV+i) (iV+f) 


16A^2 

■un f {u) f {v)dudv 
1 





r(a-l) 


0 ^0 




r(a-l) 


8A^ 




(a-1) 


2 \ , , / .,2 \ 


U- \ ^(a-1) 

8N ' 


V 

m) 


00 poo 


16A^ 



Ja-l{t)dt 


'0 JO LJO 

+ O (N-^) , N ^00. 


Ja-l(t)dt 


LJo 


Ja-l(t)dt — 1 


J L^o 


Ja-l(v)u f'{u)f{v)dudv 


The tenth term. 



uv f'{u)f'{v)dudv 

°{w)' 
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Therefore, 


^\x,y)B^^ ^\y,x)f{x)f{y)dxdy 
B^°‘~'^\x,y)Ja-i{x) \ [ Ja-i{z)dz f{x)f{y)dxdy 


f'oo poo r px 


+ -y J y Ja-i{z)dz^ y Ja-i{z)dz Ja-i{x)Ja-i{y)f{x)f{y)dxdy 

1 i*oo POO r POO px 

~ m J J J B^'"~^\y,z)dz- J B^^-^\y,z)dz x f'{x)f{y)dxdy 

1 POO POO ~ PX " " py 

~ SN J j j Jc.-i{z)dz-l j Ja-iiz)dz X f'{x)f{y)dxdy 


r*oo POO r POO 


'0 ^0 Ux 

POO POO r PX 


5^“ ^\y,z)dz- / 5^“ ^\y,z)dz / Jc,-i{z)dz Ja-i{y)x f{x)f{y)dxdy 


'0 JO uo 


Ja-i{z)dz \ Ja-i{z)dz \ Jc,-i{z)dz-l Ja-i{y)x f'{x)f{y)dxdy 


+ 0{N-^), N^oo. 

Now we want to see the mean and variance of the linear statistics ^ need 

to obtain the coefficients of A and A^, firstly we know 

/ (yWx^ ^ -XF + yF^ (ywx^, 

then we replace / with —XF + in the expression of Tr T and Tr T^, similar to previous 
discussions, denote by and the mean and variance of the linear statistics ^ {V 

we have the following theorem. 

Theorem 3.14. As N ^ oo, 


(LSE, a) _ 1 {LUE, o-l) _ 1 

dN ~ 2 4 


'0 LJO 


Ja-i{y)dy Ja-i{x)F{x)dx 


F*'" ^\x,y)dy — / B^°‘ ^\x,y)dy xF\x)dx 


Ja-i{y)dy / Jc-i{y)dy-l x F'{x)dx + O [N ‘^) 
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V 


{LSE, a) _ 


N 


1 y.{LUE,a-l) _ 1 

2^^ 4 


'0 Uo 


Ja-i{y)dy 


Ja-l{x)F‘^{x)dx 



^\x,y)Ja-i{x) 


Ja-l{z)dz 


Ja-l{z)dz 


16N 


J L-'O 

B^^-^\x,y)dy - [ B^^-^\x,y)dy 


'0 uo 



Ja-iiy)dy 

B^--^\x,y) 

B^--^\x,y) 


Ja-i{z)dz F{x)F{y)dxdy 

Ja-i{x)Jc-i{y)F{x)F{y)dxdy 
X F{x)F'{x)dx 
X F{x)F\x)dx 


J uo 


Ja-iiy)dy - 1 


B^'^-^\y,z)dz- / B^'^-^\y,z)dz 


Ja-l{z)dz — 1 


Uo 


>0 

ry 


X F\x)F{y)dxdy 


Ja-i{z)dz 


Uo 


x F\x)F{y)dxdy 


B^^-^\y,z)dz- / B^^-^\y,z)dz 


Ja-l{z)dz 


UO 


Ja—1 iy') 


X F'[x)F{y)dxdy 


1 


oo /*oo 


32iV 



Ja-l{z)dz 


>0 JO Uo 

o (Ar-2) , 


Ja-l{z)dz 


UO 


Ja-l{z)dz — 1 


J UO 


Ja-i{y)x F'{x)F{y)dxdy 


where and for N ^ oo are given in jl-fl ) and respectively with a 

replaced by a — 1. 


4 The orthogonal ensembles 

4.1 General case 

For the orthogonal ensembles, /5 = 1, the equation fll.ll) . becomes, 

r ^ 

= n \xj - Xk\Ylw{xj)[l +f{xj)]dxj, 

J laM'^ _ _ -AT 1 


l'^j<k<N 


i=i 


and we assume N is even. Here, 

41) _ 


_ 

Cjv — 


f(a,b)^ ni<j<fc<Ar l^i ^k\ 11^=1 w{Xj)dXj 


depends on N. 
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We also follow the treatment of [6l |T5], firstly using Theorem 12.51 and some computations, we 


hnd 




= det 


b pb 



£{x,y)T^j{x)Tik{y)w{x)w{y){l + f{x)){l + f{y))dxdy 


N-l 


j,k=0 


where is a constant depending on N and 'n'j{x) is an arbitrary polynomial of degree j. Let 


'lpj{x) = 7lj{x)w{x), 


(4.1) 


it follows that 


GSi^(/)]' = det(/+(M«) 'lW), 


where 


= ( [ 'ipj{x)e'ilJk{x)dx'\ , 
a / j,k=0 

= ( / {f{x)ipj{x)eiJk{x) - f{x)'ipk{x)ei;j{x) - f{x)'ilJk{x)e{f'ipj){x)) dx 


N-l 


j,k=0 


If ^ = {yjk)fJo^ then 


where is an integral operator 


G^i\f) =det{I + K\^>f), 


( 1 ) 


rW _ 


- Eyfc=o djki^j ® ^i’k T.j,k=o djki^j ® i^k 
- E^a^o djk^i^j ^e'tpk-e EyfcEo djkS^j ® i’k^ 


^d)2) 
iy(2>l) ^(2,2) 


In the next theorem, we obtain relations on 


Theorem 4.1. 


I^(l-2) _ ^(2,1) _ /^(2.2) _ 

^N,l ~ 1 J^N,1 ~ X^N,1 ^ ^N,l ~ t. 
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Proof. For any integrable function g[x) supported on [a, 6], we have 

DK'Sfg(x) = P f KSf(x,y)g(y)dy 



= J -^K^N,i\^^y)9{y)dy 

P Q 

= ^'^ji^)9jk'^k{y)9{y)dy 

^ j,k=0 

l>b 

= / Y^d9e'ijjj{x))gjk'tpk{y)9{y)dy 

d°- j,k=o 

pb N-1 

= Y '^j(^)Pjkf^k{y)g{y)dy 

j,k=0 

= [ K^i?{x,y)9{y)dy 

J a 

which implies, 

^(l>2) _ ^^(2,2) 

^N,l ~ ^^N,l ■ 

Note that 

N-l 

= Yj yjksf’j ® i^k£ 

j,k=0 

N-l 

= ~Y ^^i^k (4.2) 

j,k=0 

that is, 

iv(2.1) _ /v(2>2) _ 

Moreover, from fj4.2p. 

N-l 

DK^^fe = - gjk{Dei>j) 0 eipk 

j,k=0 

N-l 

= - ^ gjk'tpj (g) eifk 

j,k=0 

— ^N,l ■ 

The proof is complete. 

□ 
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The series of computations presented below takes the determinant into the form of identity plus 
scalar operators. From Theorem 14.11 can be written as 


DK'^fs 



DK 


( 2 , 2 )' 


Af,l 

AtF - . 



where 


A := 


By Theorem 12.61 we have 


D 

o\ ~ / 



VO 

I \ 


^N,l ^ 


K 


e K 


( 2 , 2 ) 

N,1 

( 2 , 2 ) 

N,1 


(/) = det ( / + ABfJ = det (^/ + BfAj . 


Since 


BfA = 




K 


( 2 , 2 ) 


N,1 


iv(2.2) ^(2,2) 




K 


N,1 




e K 


( 2 , 2 ) 


N,1 



K'^^fefD A'lf/ 
KSfefD - efD 


g2>(/) 


= det 


then 

I + K^^fefD K^^ff 
K^^fefD-efD I + f^ 

The computation below reduces the above into a determinant of scalar operators. We subtract 
row 1 from row 2, 

G'Af) 


1 2 


= det 


I + K^j^fsfD 


-I - efD 


I 
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Next, add column 2 times / + efD to column 1, 



/ + KSfefD + + efD) f 

0 I 

det (/ + K^^fiefD + f + fefD)^ . 


,( 2 , 2 ) 


,( 2 , 2 ) 


det 


Remark. The above result agrees with [6] for the GOE case if we take / = —/r xj^ where xj is the 
characteristic function of the interval J. 


Using fl3.3p . Lemma 12.71 and note that / is a Schwartz function, then we have the following 
theorem. 


Theorem 4.2. 


g!)’(/)]" = det (/ + A')):? (f + 2/) - A<Jf £/' - A'JJ.i'/e/') . 

where the kernel of reads 

N-l 

^ tijke'fj{x)f)k{y). 


j,k=0 


4.2 GOE 


It is convenient in this case to choose 

/ \ — — TTTi 

w{x) = e 2 ^ X G M, 

following Dieng and Tracy-Widom’s discussion P[T5]. We want to choose 'ipj to make simplest 
possible. Dehne 

d 

'<p2n+l{x) -.= —ip2n{x), 1p2n{x) := ip2n{x), U = 0, 1, 2, . . . , (4.3) 

ax 

where ffj{x) is given by (I3.5p . 

Cj = 7r3 22^r(j + 1). 

3.2 

We can check that this dehnition satishes fl4.ip . i.e., fjj{x) = 7rj{x)e~^, j = 0,1,2,..., where 
'Xj{x) is a polynomial of degree j. The matrix ;= ( 7pj(x)e'f>k(x)dx) is computed below 

\ / j,k=0 

from fl4.3p . 
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Lemma 4.3. (Dieng, Tracy — Widom) 

0 1 0 0 


= 


-10 0 0 
0 0 0 1 
0 0-10 


0 o\ 
0 0 
0 0 
0 0 


0 0 0 0 ■■■ 01 

\o 0 0 o--- -loy 

It’s obvious that ^ = —M^^\ so H 2 j, 2 j+i = —l 5 h 2 j+i, 2 j = 1, and fijk = 0 for other cases, 


NxN 


hence 


N-l 




j,k=0 


—-1 


—-1 


Sij2j{x)ij2j+l{y) + Y ^^2j+l{x)ij2j{y) 


j=0 j=0 

N 1 N 


i\ _1 i\_1 

2 ^ 2 ^ 

Y '^2j{x)<^2j{,y) - Y 

j=0 j=0 


By fIXTTD . we have 


iL-i 

n X 


—-1 


Y ^^2jix)ip2j{y) = Y ^/jy^2j-iiy) - \ J + 7:^2j+iiy) 


j=0 


j=0 

^-1 


41-1 


Y Vj^^2j{x)(p2j-i{y) - \/i + 7^£y^2jix)^2j+iiy) 


j=0 

^-1 


j=0 

N 


Y Vj^y^2j{x)(p2j-l{y) -Y\H 7^£^2j-2ix)^2j-l{y) 

3=^ 


3 = ^ 

N 


^ ^ JSJ' 

Y Vj^V2j{x)^2j-l{y) - Y —e^N{x)>^N-i{y) -YY~ 

i=i ^ i=i ^ 


7 ^^‘f 2 j- 2 (x) - x/jey> 2 j(x) j y> 2 j-i(y) 
i=i 


e(fijv(x)(f>jv-i(y). 


Using fl3.11|] again, we hnd 


^2j-l(x) = \jj- ]^ip2j-2{x) - V7<d2i(x), 
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then according to Lemma 12.71 


^2j-i{x) = - ^e(p2j-2{x) - ^/jeip2j{x). 


Hence 


It follows that 


K 

2 


eip2j{x)ip'2j{y) = - ^(/92i-l(x)(/72i-l(2/) - y —eipN{x)ipN-i{y)- 


j=0 


j=l 


E. 

2 


= Y 7’2i(a^)7’2i(|/) + ^(^2i-i(a:)(/?2i-i(2/) + A/ —e^N{x)^N-i{y) 
j=0 j=l ^ 

^-1 

= Y + Y -j£¥’N{x)^N-i{y) 


j=0 


= SN{x,y) + \l—eipN{x)(pN-i{y), 


where 


e ^ ^ ^ y^N{x)ipN-i{y) - ^N{yWN-i{x) 

SN{x,y) ;= 2 ^ ^j{x)^j{y) = \J— ---. 

i=o ^ 

Here we have nsed the Christoffel-Darbonx formula in the last equality. Note that 

SN{x,y) = K^^\x,y). 

By Theorem 14.21 we have the following theorem. 

Theorem 4.4. 


G«(/) = det J + + 2/) - SMsf - S^fef + \ -s^pn ® + 2/) 


+ ® eipN-i)f + \ {ep>N ® e{p>N-if)) /'^ • 


4.3 Large N behavior of the GOE moment generating function 


Now we consider the scaling limit of G^^\f) 


. Write 


G^a) =:det(/ + r), 


where 


T-. = SN{f + 2f)-Sr,ef-SNfef+^e^M®^N-i{f + 2f) 

[n fW 

+ Y yo e(pN-i)f + Y y ® f ■ 
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In the computations below, we replace f{x) by / ( \/2Nx) and f'{x) by 


d 


f ( Vmx) = ^ / ( Vm: 


\/2Ndx 

Similarly to the GSE case, we have the following theorems. 


IX ] . 


Theorem 4.5. 


1 / X y \ ~ y) 

v™ vW VW; vr(x-|/)’ 


lim -^Sn ‘ ^ 


a; \ 1 


VmJ vr 


Theorem 4.6. 


^Pn-1 f , j = —(—l)2247r '^N 4sinx + 0 4 j ^ ]\r qq, 


VvW/ 


( X \ ,,iVll 3 5\ 

--j== \ = (—1) 2 2 47r 2 iV 4 sinx + O yN 4 j ^ jy 




=-2-2A^-2+ 0 , A^^CX), 


e ((pw-i/) 


X 


\ / X 5 1,^3 

- I = —( — 1) 2 2 47r 2 N 4 


sin?//(|/)d2/- / sin?//(?/)d?/ 


J 

Using Theorem 14.51 and Theorem 14.61 to compute Tr T and Tr as A^ —)■ 00 . we have 


+0 (^■^) , N ^ 00 . 


1 

Tr T = — [/“^{x) + 2f{x)~\ dx 

'TC } — 


27r\/^ J- 
1 




00 3; 


X - ?/ 


'—00 ^ y 


f'{x)dx 


2tiN 


sin^ X [/^(x) + 2/(x)] dx — ^ 


sinx f'{x)dx + O [N 2^ N ^ 00 , 


and 


1 /.OO /.oo 

sin(x — y) 

TrT^ = — / / 

TT^ / / 

J —00 v/ —00 

[ x-y \ 


1 2 


f*oo roo 


7T- 


sin(x — y) 
'y/m J-ooJ-00 x-y 


2y/^N J-0 

[f{x) + 2/(x)] [f{y) + 2f{y)] dxdy 
Si(a; - y)f{x) [f{y) + 2f{y)] dxdy 


poo roo 


7T- 


'y/m J- 


00 J —00 


sin(?/ — z) 
y- z 


fiz)dz 


sin(?/ — z) 


y- z 


f{z)dz 


sin(x — y) 
x-y 


f{x) [f{y) + 2f{y)] dxdy + O {N , N ^ 00 . 
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Now we want to find the mean and variance of the linear statistics ^ ( \/2A^Xj j, since 




X 


\ 2 

-XF [V2Nx^ + , 


we replace / with —XF + in the expression of Tr T and Tr T^, we have 


log det (/ + T) 

2 ^oo ^ noo 

= —Xi — / F{x)dx -— / sin^ a; F(a;)dx 

Tr —r^ ttjV 


- 1 ) 




^(4 
2 I^Tr 
1 


oo ^ 

F‘^(x)dx - - 


— OO 
poo 


poo 

sin(x — y) 

J —oo 

x-y 


TT 


VM J-c 


r*00 /‘OO 


sin(x — j/) 
x-y 

sin(x — y) 


F{y)dy 


2y/^N J- 
1 2 

F{x)F{y)dxdy 
sin(x — y) 


sinx F\x)dx + O [N 2 


x-y 


-F{y)dy 


F\x)dx 


'^i^x — y)F\x)F{y)dxdy + O (yN >, N ^ 00 . 


J-c 

Denote by and the mean and variance of the linear statistics ^ ^\/2A^a;j j, then 

we have the following theorem. 

Theorem 4.7. As N ^ 00 , 


(GOE) _ (GUE) 

En ~ rw 


2eN 


sin^ X F{x)dx 


-ly 


4:\/^N J-c 


V 


(GOE) 

N 


= 2V 


(GUE) 

N 


2 nV^ J-00 VJx x-y 

sin(a: — y) 


foo poo 


-00 x-y 

Si(x — y)F'{x)F{y)dxdy + O {N~y 


sinX F'{x)dx + O [N 

F\x)dx 




J- 

where and for N ^ 00 are given in M.X\) and M.3^) respectively 


4.4 LOE 

We now specialize results obtained in Section 4.1 to the situation where w is taken to be the square 
root of the Laguerre weight, namely, 


q : X 


w{x) = x‘ 2 e 2 , a > —2, x G M"^. 

Again we choose a special tjjj so that is as simple as possible. Let 

i^ 2 n+i{x) := ^p 2 n{x) := u = 0, 1, 2, . . . , 

dx 


(4.4) 
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where and are given by 


(oi+l) / \ 

if) \x) := 




a;2+ie 2, j = 0,1,2, 


d«+i) 


(4.5) 


It’s easy to see that 


T (t\ 

:= xte~t, j = 0,1,2,.... 

S’ 


= ^fc, j. A; = 0,1,2, 


The next theorem shows that fl4.4p satishes fl4.1l) . namely, 'il)j{x)x~^e^ is a polynomial of degree j. 
Theorem 4.8. 'ijjj{x)x~^e^, j = 0,1,2,... is a polynomial of degree j. 

Proof. We prove this by considering two cases, j even and j odd. It is clear for even j, that 
'f’ 2 j e®/^ is np to a constant mnltiple of (x). 

If j is odd, then we hnd. 


ij2n+l{x)x 2 62 = 




_ ^ ^ 

X 2 62 


d“+l) 


r(a+l)/ \ “4-1 -£ \ £ 

L^n [x)x‘^^ e 2 X 2 62 


^2n 


^2n 


X ( L 


ir'l-t))' + (f +1) iT'h) - 


, n = 0,l,2,.... 


Thus 'tjj 2 n+i{x)x~°'^‘^ eS is a polynomial of degree 2n + 1. The proof is complete. □ 

We use fl4.4p to compute := fjj{x)e'ipk{x)dx'j^f^, resulting in the following theorem. 

Theorem 4.9. 


0 1 0 0 
-10 0 0 


Mi^'> = 


0 0^ 

0 0 


0 0 0 1 ■■■ 00 

0 0 -1 0 ■■■ 0 0 


0 0 0 0 
y 0 0 0 0 


0 1 
-1 0 / 


NxN 
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Proof. Let rrijk be the (j,/c)-entry of again separating into four cases: {j,k) = (even, odd), 

(odd,even), (even, even) and (odd, odd). We find, for the (even, odd) case, 


'nr2j,2k+i 


'i(j2j{x)e'ilj2k+i{x)dx 




~(a+l) 

P2j 


x)pi'^^\x)dx 


= S. 


jk- 


For the (odd, even) case, note that is antisymmetric, then 

^2j+l,2k = —'m2k,2j+l 

— ~djk. 

The (even, even) case, 


Xn2j,2k — 


f)2j{x)ef)2k{x)dx 


'0 


„ ,...| ,...| / 4r’WxSe-i 

2cL. TL do wo 


-'2j ^2k 

1 


2d?/- / (?/)?/2e ) dx 


2cr‘'cr‘'^o '' 



4“+^^(x)x2e 2 (2 / L^“+^^(?/)?/2e / L^"+^^(l/)l/2e 2(iyMa: 


9„(“+L («+i) 
^2fc 



L2“’'’^^(x)x2e 2dx L^2t^^\y)y'^^ 


L^^'^^\x)x^e~^dx / L')2]^^\y)y^G~^dy 


2“+2r (j +1 +1) r (/c +1 + f) 2“+2r (j +1 +1) r (A: +1 + f) 


9„(“+b (a+1) 
^^2i ^2fc 




j!A;! 


= 0 , 


where we have used the fact 


L2“^^^(x)x2e 2dx = 


22 


0 


, J = 0,l,2, 
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Finally, the (odd, odd) case. If j < k, 

poo 

m2j+i,2k+i = / ip2j+i{x)e'ip2k+i{x)dx 

Jo 

r 1 ' 

= eDip^^^'^^\x)dx 

r 1 ' 

= / (f^^^'^^\x)dx 


10 a 


(a+l) 


2i 


^ \^)) + ( ^ + 1 j -^2j 


a 


(a+l): 




^ 

X 2 e 2 


2cix 


O ^ O y 


-2k 


(a+l) (a+l) 

^2j ^2k ‘J 


4d‘F) 


I (Fy’w) + (f+1 


L. 


r‘’w-2='dr’w 


a;“+^e“ 


= 0 , 


since, 2j + 1, the degree of the polynomials x + (f + l) ~ is less 

than 2k. 

If j > k, due to the fact that is antisymmetric. 


^2j+l,2A:+l — ~^2A:+l,2j+l 

= 0 , 


hence 


Tn 2 j+l, 2 k+l — 0, j, /c — 0, 1, 

It is the desired result for The proof is complete. 


□ 


We begin here a series of computations analogues to those in derivation of the LSE problem, 

l2 


ending up with an expression for 

, -1 




as a scalar Fredholm determinant. 

It’s obvious that ^ = —M^^\ so ^ 2 j, 2 j+i = —l,h 2 j+i, 2 j = 1, and fXjk = 0 for other cases, 

hence 


Af-l 




j,k=0 


—-1 
2 


—-1 

2 


- S^/J2jix)^/J2j+liy) + Y ^'^2j+lix)'llj2jiy) 


4-1 


Y - Y [p2Y\y) 

j=0 j=0 


^dx 
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Differentiating fl4.5p with respect to x, we find, 




d“+i) 


x(L 




^ 

x2e 2. 


Recall that the Lagnerre polynomials satisfy the differentiation formnlas [8] 

- {j+ a + l)L^p^^\x), j = 0,1,2,..., 

X = {j + - (i + a + 2 - x)L^^^^\x), j = 0,1,2,.... 

The snm of (I4.6p and (14.7p divided by 2, gives, 

^ (iT'T))' = 3 = 0.1,: 

which is the same as, 

X (Lf«'(x))' + Lfy‘>(x) - Lf»\.). 

Hence the derivative of !.p^^^^\x) becomes, 

(o+l) / \1 
Lp) ’[X) = 


d«+l) 




a X 

X 2 e~ 2 


^V(i + l)(i + « + 2) + a + l) (p5“t^^(x). 


So replacing j by 2j, we see that. 


2J y 2 
Continning, we compnte the snmmation to give. 


^-+H (J+-1 d 2+^) dtf (!/)■ 


(4.6) 

(4.7) 


(4.8) 


j=0 


f-i 


El 


j-o 





r(«+l). 


X 
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From fl4.8p . and using Lemma [2.71 we have 





~(q:+1) 


x) 




The sum simplihes to 


[fiT^yy) 

j=0 


E. 

2 


+ -VN{N + a + l) 


i=i 


The computations above gives a compact form for y) : 


KN’?i^,y) = Y1 + + a + 1) 


E 

2 


("+1)/'^',-(«+!)/ 


-(a+l) / \ ~(a+l), 


j=0 

N-l 


i=i 


Y1 - -\/N{N + a + l) eip^^^^\x)ip^^^l\y) 


j=0 


= SN{x,y) - -VN{N + a + 1) eip^^^^\x)^^^^l\y), 


where 


N-l 

SN{x,y) := ifY^\x)^Y^\y) = -^/N{N + a + l) 

j=0 




x-y 


Here we have used the Christoffel-Darboux formula in the last equality. 
By Theorem 14.21 we have the following theorem. 


Theorem 4.10 

41) 




= det + SNif + 2/) - S^ef - SNfef - -^/N{N + a + l)[eyp^Y'^ 0 IpYi {f + V) 
- ]^^N{N + a + l) 0 e^Yi) f - \^/N{N + a + l) 0 e (^KV)] f 


4.5 Large N behavior of the LOE moment generating function 


Now we consider the scaling limit of 


g"(/) 


We write 


GiJ’a) =:det(/ + r), 
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where 


T: = S^if + 2f) - S^ef - S^fef --^N{N + a + 1) ® {f + 2f) 


-^N{N + a + l) ® f 


- -^N{N + a + l) 


/'. 


In the computations below, we replace f{x) by / ANxj and f'{x) by 

d 




\/ANd-^/x'‘ 

Similarly to the LSE case, we have the following two theorems. 


Theorem 4.11. 


lim — 

V—>-00 2,N 




N^oo \AN' AN 


x2 J.2 


where {x,y) and B^°‘~^^^{x,x) are given by U.di) and (fi.?|j with a replaced by a + I- 


Theorem 4.12. 


x'^ \ 




i '^q;+i(^) 


X 


+ 0(^N-l^ , N ^oo, 


'(a+l) / ^ 

(iiV 


= N- 


■Ja+iiy)dy -1 


L ./0 


+ 0 ( 


N' 


N —>■ oo. 


~(q+1) 

^^N-1 


/ X^ 

V 4iV 


= N- 


J^+i{y)dy + 0(N 


N 


CX), 


£ 




2 


Ja+i{y)f{y)dy - / Ja+i{y)f{y)dy 


+ 0 N 


N 


oo. 
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Using Theorem 14.111 and Theorem 14.121 to compute Tr T and TrT^ as N ^ oo. We obtain 


TrT = 


^ ^0 
1 

1 

4iV 


I [f‘^{x) + 2/(x)] dx 

0 

1 

I Ja+i{y)dy-l 
0 

oo nx 

B^'^^^\x,y)dy — / B'^°'^^\x^y)dy 
0 IJx Jo 

OO nx 

B^^^^\x,y)f{y)dy- / B^°^^\x,y)f{y)dy 

Jo 


Ja+l{x) [f{x) + 2f{x)] dx 

X f'{x)dx 


>0 

1 r°° 

1 r°° 

m 


X f'{x)dx 


Ja+i{y)dy - 1 


Ja+i{y)dy - 1 


Ja+i{y)dy 


'0 Uo 
+ O (N-^) , N ^oo, 


L-'O 


L-'o 


X f'{x)dx 


Ja+i{y)f{y)dy - / Ja+i{y)f{y)dy 


X f{x)dx 
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and 


[f{x) + 2f{x)\ [f{y) + 2f{y)\ dxdy 

S("+^)(x,|/)J«+i(x) [ Ja,+i{z)dz-l [f{x) + 2f{x)][f{y) + 2f{y)]dxdy 

Jo J 


^CXD poo r px 


'0 Jo L^o 


Ja-\-l(^^^dz 1 / J(yj^i(^Z^dz 1 t/a+l (^) t/a-|-l (y) 


[f{x) + 2f{x)] [f{y) + 2f{y)] dxdy 

1 POO POO r POO PX 

— / B^'^^^\x,y)\ B^'^^^\y,z)dz- B^^^^\y, z)dz x f{x) [f{y) + 2f{y)] dxdy 


— I B^^^^\x,y)y B^^^^\y,z)f{z)dz-B(^^^\y,z)f{z)dz 

X f{x) [f{y) + 2f{y)] dxdy 

1 POO POO r PX "I r py 

— / / Ja+iiz)dz \ Ja+iiz)dz-l X f (x) [f (y) + 2f (y)] dxdy 


Zn j j [y Ja+i{z)f{z)dz - j Ja+i{z)f{z)dz^y Je,+i{z)dz-l 

X fix) [fiy) + 2 fiy)] dxdy 

1 POO POO r POO PX “I r px 

J J “ y B^'^^^\y,z)dz J Ja+lifdz -1 Ja+fy) 

X fix) [fiy) + 2fiy)] dxdy 

1 POO POO r POO PX “I r PX 

— \ B^'^^^\y,z)fiz)dz- B^‘^^^\y,z)fiz)dz \ Ja+iiz)dz - 1 Ja+fy) 


'0 Jo Ux 


X fix) [fiy) + 2 fiy)] dxdy 


r>00 POO r PX 


'0 JO UO 


Ja+liz)dz-l / Ja+liz)dz-l / Ja+liz)dz Ja+liv) 


X fix) [fiy) + 2 fiy)] dxdy 


"OO POO r rx 


'0 JO IJO 


Ja+liz)dz-l / Ja+iiz)dz-l 


Ja+iiz)fiz)dz - Ja+iiz)fiz)dz Ja+iiy)x f ix) [f iy) + 2f iy)] dxdy 


+ 0{N-^), N^oo. 

Now we want to see the mean and variance of the linear statistics Ylf=i 
to obtain the coefficients of A and A^, hrstly we know 

/ (yw^'^ ^ -XF (yify + yF^ (yzfy, 


Jxj ), so we need 


then we replace / with —XF + in the expression of Tr T and Tr T^, similar to previous 



discussions, denote by and the mean and variance of the linear statistics Y^^=i ^ (a/4A^Xj) , 

then we have the following theorem. 


Theorem 4.13. As N ^ oo, 


{LOE,a) 

EN 


{LUE, a+1) 


1 f°° 

^ Jo 

1 

m ^ 


Ja+i{y)dy - 1 


^ ./o IJO 
B^'^^^\x,y)dy- [ B^^^^\x,y)dy 


Ja+l{x)F{x)dx 

X F\x)dx 


0 Uo 


Ja+i{y)dy - 1 


'0 

rx 


Jo,+i{y)dy 


J Uo 


X F'{x)dx + O [N 


,,(LOE, a) 


c\ ^^{EUE^ a+1) 


>0 


Ja+i{y)dy -1 
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where and j\r qq given in ljil.4\ ) and U.5\) respectively with a 

replaced by a + 1. 
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5 Conclusion 


In this paper, we study the moment generating function of linear statistics of the form ^{^ 3)1 
namely, the expectation value or average of 

N 

exp -A^F(a;j) 
i=i 

computed in the “background” of the symplectic and orthogonal ensembles. We then specialize to 
the Gaussian case, where the background weight is the normal distribution and the Laguerre case, 
where the background weight is the gamma distribution. 

Finally, we compute the large N behavior of the moment generating, where F{.) is suitably 
scaled and obtained the mean and variance of the corresponding linear statistics in the four cases. 
The more complex situation with the Jacobi background will be studied in the future, in which case 
the polynomials depend on two parameters a and b. 
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